Stability and correlations in dilute two-dimensional boson systems 
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The hyperspherical adiabatic expansion method is used to describe correlations in a symmetric 
boson system rigorously confined to two spatial dimensions. The hyperangular eigenvalue equation 
turns out to be almost independent of the hyperradius, whereas the solutions are strongly varying 
with the strength of the attractive two-body potentials. Instability is encountered in hyperangular, 
hyperradial, and mean-field equations for almost identical strengths inversely proportional to the 
particle number. The derived conditions for stability are similar to mean-field conditions and closely 
related to the possible occurrence of the Thomas and Efimov effects. Renormalization in mean-field 
calculations for two spatial dimensions is probably not needed. 
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Introduction, howei dimensions than three are nec- 
essary in several branches of physics, for instance sur- 
face physics, semiconductor physics, artificial atoms, and 
quantum dots. The advanced tools of well-controlled ex- 
ternal fields, used in atomic and molecular physics to 
manipulate the effective interactionsLallow confinement 
of the systems to lower dimensions 1^. The experi- 
mental investigations employ continuous variation of the 
dimension by use of tunable deformed external fields. 

The basic properties vary dramatically with the dimen- 
sionality of the system, as highlighted in two dimensions 
(2_D) where the centrifugal s-wave barrier is negative for 
two particles and even an infinitesimally small attraction 
provides a bound state 0, ■ The stability is strongly 
dependent on the deformation or the effective dimension 
of the confining potential 0, 0, ■ The simplest A^-body 
structures are the Bose-Einstein condensates for identi- 
cal bosons. These systems are dilute, weakly interacting, 
and well described by mean- field models '3,^ with an in- 
teraction strength adjusted to reproduce the low-energy 
scattering cross section in the Born approximation. The 
interaction in 2D is then obtained from short- or zero- 
range three-dimensional (iP) potentials restricted to the 
mean- field Hilbert space |^|, |0| . 

Inclusion of correlations in the wave function prohibits 
this renormalization. Instead a finite short-range poten- 
tial with the correct scattering length should be used. 
If the correlations are appropriately accounted for, the 
large-distance behavior must come out correctly with the 
realistic interaction The consequences for lower 

dimensions are not yet investigated. The huge difference 
between two- and three-body properties in two and three 
dimensions is most likely more pronounced for A^-body 
systems. Experimental results from varying dimension- 
ality are easier to interpret if the limits are known. 



The purpose of this letter is to (i) formulate a frame- 
work for investigations of correlated structures under 
strict 2D confinement for N interacting identical bosons 
in an external harmonic field, (ii) derive stability con- 
ditions in terms of particle number and two-body inter- 
action properties, (iii) extract the basic features of the 
solutions, and (iv) provide a conceptual link between the 
successfully renormalized mean-field models and the cor- 
related solutions with the bare (effective) interaction. 

Theoretical method. We shall briefly sketch the 
method to establish the notation and the pertinent for- 
mulae for two dimensions. We follow the derivation for 
three dimensions given in The A^ identical bosons 

have masses m and coordinates r^. We use the hyper- 
spherical adiabatic expansion method where the only 
length coordinate is the hyperradius p defined by 



1 ^ 
- Vr2. 

i<j 



N 



(1) 



where fij — ri 
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and R — J2i ^i/^ is the center-of-mass 
coordinate. The remaining 2N — 3 relative coordinates 
in two dimensions are angles where we define cty related 
to the size of by = V2p sin . If permitted in the 
context, we shall omit the indices ij. 

The center of mass separates out and we only need 
to deal with relative coordinates. The related volume 
element is p^^~'^dpsinacos^^~^ adadi9dilAr_2, where 
a = aij, d describes the direction of and O.n-2 de- 
notes the remaining angular part of the volume element 
corresponding to the last N — 2 relative vector coordi- 
nates. 

An external harmonic potential rriLo^ rf /2 of angu- 
lar frequency lu is by use of eq. divided into a center- 
of-mass part and a hyperadial part. The relative Hamil- 
tonian is then separated into a hyperradial part and a 



2 



hyperangular part, ha, i.e 

fi2 / 1 



H = 



2 m 



^ 2N-3 ^ 



dp 



P 



dp bf ' p 



2 



(2) 



where G = G{T,a) is linear in both and V and 
the reduced potentials are given by 



(2iV-5)(2iV^7)^^^,^_ 



cot^ a 



47V- 9 



(12) 



angle 



o 2 N 

i<j 



2N-6 - {2N - 4) cos 2a d 

dcP- 



sin 2a 



dc 



(3) 



(4) 



where all differential angular dependence, except a, is 
collected in Dangle- The trap length ht is given by 
6( = h/(muj). The two-body interaction V^j is of short 
range, e.g. a Gaussian Vq exp(— r^^ /6^) or a square well 
Vb6(?'ij < 5), where is the truth function. 

The relative wave function ^(p, ri) obeys the 
Schrodinger equation 



H'9{p,n) ^ E'^ip,n) , 



(5) 



where E is the energy. We write as an adiabatic ex- 
pansion 1^ 01 where the first term is 



^ip,n)^p-^^^~'y'f{pMp,n) 



(6) 



with the hyperradial volume element explicitly extracted. 
The angular wave function ^{p, fl) is for fixed p an eigen- 
function of hfi with the eigenvalue X{p), i.e. 



hn<fip,n) ^ x{p)^p,n) . 

The corresponding radial equation is then 



(7) 



(8) 



2mC/(p) _ A(p) (2jV-3)(27V-5) 



where the adiabatic potential C/ is a function of the hy- 
perradius consisting of three terms, i.e. the angular aver- 
age A of the interactions and kinetic energies, the gener- 
alized centrifugal barrier, and the external field. 

For large particle distances only relative s waves con- 
tribute. With a Faddeev decomposition of the angular 
wave function only the dependence on distance a.y is 
left, i.e. 



N 



N 



Hp^ aij) 



i<j 



- - M— H 

j^j- sm2 aij cos 2 



(10) 



where we again explicitly extracted the square root of the 
volume element. 

The integro-differential equation for (jj^p, a) is obtained 
from eq. (jT)) by integrating over all other angles than 
ai2 — a, denoted by r TiJ, i.e. 



- + Vc + V + vi + V2 - 4>{p, a) = y '^'^ ^ '^^^ 



v{a) = 4aB(^)^ exp(-2(psina/6)^) , 
v,{a) = 2aB{^)\N-if{N-2) 

X J dxx'^-*exp[2ix-l)i^co8af] (14) 

. i^~^nN~2)a, 
cos^ a 

V2{a) « ^wi(a)(l - itan2a)^~''e(a < 7r/3) , (16) 



(13) 



as 



' V{r)d'drdr ^ 



2'nK^ 



2h'^ 



(17) 



The last expression for as is valid both for Gaussian and 
square-well potentials V. The approximations for vi and 
V2 are very accurate for pcosa » b. For a — n/2 we get 
exactly vi(7r/2) = 2{N - 3){N ~ 2)aB{plbf, and V2 is 
very small. 

For a short-range interaction the right-hand side of 
eq. (|ll|l is independent of p as well as vi and V2 when 
a is not too close to tt/2. The only p dependence is then 
through V, which approaches a zero-range interaction in a 
as p increases. The eigenvalue \{p) is therefore expected 
to be constant in large ranges of p. These features are 
unique for two dimensions. 

Stability conditions. For comparison we first consider 
the mean-field approximation in two dimensions. This 
was investigated in details in |l5j | with a potential de- 
rived from a three-dimensional zero-range potential by 
lla. N eglecting the logarithmic energy dependent term 
in |l5l Ha l only a two-dimensional zero-range potential 
remains precisely as the delta-function limit of our Gaus- 
sian short-range potential. In general, for an interaction 
of short range, i.e. small 6, the differential equation is in 
dimensionless quantities given by 



ld_ 

X dx 



2as(^-l)|/m(a;)P-e /™(a;) = 



(18) 



where x = rjb and e = 2mEmb^ / h? are measures of the 
single-particle mean-field coordinate r and energy Em- 
The radial wave function is approximated by a Gaus- 
sian, i.e. fra = ex.p[—x'^/(2d?)]/{y/7rd) renormalized as 
/ d-&xdx\fm{x)\'^ = 1 as in [13. The corresponding en- 
ergy per particle e is then as a function of d given by 



d^ 



1 



b 

which only has a minimum when 



-[l + aeiN-l)] 



aB{N-l) > -1 . 



(19) 



(20) 
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Then the energy and the Gaussian width are 



b 



[1 + (A^ - l)aB 
6 \2 



1/4 



>(-) ^l + iN~l)aB. 



(21) 
(22) 



These results coincide for TV — 1 w and as = g with 
those derived in If the interaction strength in an 

experiment suddenly is changed from repulsive to an at- 
tractive value as, the new state of the A^-body system 
can only be stable for particle numbers smaller than the 
critical value Nc = 1 — 1/as- For N > Nc the motion 
is towards larger densities and either a total collapse or 
a reduction of particles in the gas caused by molecular 
recombination. 

Thus, only repulsive or very weakly attractive poten- 
tials provide stable mean-field solutions. The basic rea- 
son is that even infinitesimally small attractions bind two 
particles in two dimensions. The corresponding two-body 
Schrodinger equation with wave function /2 and energy 
E2 is 



ld_ 

X dx 



2aB 



V{x) - E2 



Vo 



f2{x) = , (23) 



which for a weakly attractive interaction has the bound- 
state energy E2 = —Ah?/{m.lP')exp{2/aB — 27) and the 
mean-square radius (r^) ~ 2fi,^/(3m|i?2|), where 7 is Eu- 
ler's constant, see 
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FIG. 1: The centrifugal terms Vc and vi from eqs. (1121 and 
IIH as function of a. The thick lines show Vc + vi for TV = 20, 
various values of as, and p/b = 10'' unless otherwise indi- 
cated. The thin, solid lines show the contributions from vi 
and Vc for as ~ —0.02. 

The radial potential in eq. @ depends crucially on A 
determined from eq. I|ll|) which in turn is dominated by 



the terms Vc and wi shown in fig. ^ When a approaches 
7r/2 the approximation in eq. I|15|l and Vc eq. (I12|) both 
diverge as cos~^ a « {t^/2 — a)^^. Thus, when qb < 
~{N - 5/2){N - 7/2)/[{N - 2){N - 3)^] « -1/(A^ - 2) 
an attractive pocket inevitably appears for large p. This 
divergence at the point a = tt/2 disappears when the 
exact expression in eq. H14|l is used for vi . Then the low- 
est eigenvalue A would be finite and proportional to p~^. 
These solutions correspond to many particles close to- 
gether which probably violates our assumption of s-wave 
dominance in the wave function. To avoid this divergence 
at large p the strength of the attraction must be limited 
by aB{N —2) > — 1, which is almost identical to eq. H20|l . 
Thus, remarkably enough the mean-field stability condi- 
tion is precisely also obtained from the angular potential 
when N -2 7ii N. 

For small a the angular centrifugal term Vc from 
eq. (|12|l diverges as — a^/4. This is the limit rigor- 
ously separating attractions leading to either no bound 
states or infinitly many bound states of Thomas or Efi- 
mov character 3]. Thus, a small two-body attraction v 
is for large p sufficient to bind a state in the pocket at 
small a. This is seen by substituting x = \/2pa/b in 
eq. ()11|1 when only the potential v and Wc are included. 
The energy and the mean-square radius of the solution 
becomes A(p) 2mp^E2/T? and {o?) 2/(3|A|) = 
l/12(&/p)^ exp(27 — 2/05). If the size in a space has 
to be smaller than unity, p/b must exceed exp(— l/a^) 
which is huge when —1/N < ub < 0, i.e. the interaction 
is attractive but allows physical solutions in agreement 
with eq. Thus, the diverging A corresponding to 

the bound two-body state is never encountered because 
either the interaction is too attractive leading to solu- 
tions in the pocket at large a, or the interaction is too 
weak to bind at small a for p values less than the trap 
length. Therefore diatomic recombination is unlikely. 

In any case the angular potential provides two types 
of minima at small and large a-values, respectively. The 
large a-behavior corresponds to the stability condition 
obtained by mean-field calculations. The structure at- 
tempts to maximize the two-body attractions by a rather 
similar distance between all particles confined to the vol- 
ume allowed by the given hyperradius, i.e. two particles 
are far apart but all others are correspondingly close. 

The small a-behavior turns out to be unimportant for 
stability, because the extremely weak binding requires 
a spatially extended wave function attempting to reach 
beyond the confining boundaries of the trap. A huge 
trap would in principle allow binding by this pocket, 
which corresponds to a structure with one two-body 
bound state and consequently this does not resemble a 
condensed state. Such a structure would not be found 
by mean-field calculations although the two-body bound 
state always is present for overall attraction in two di- 
mensions. 

The deeper-lying reason for not populating this two- 
body bound state seems to be related to the choice of 
strength for a two-dimensional zero-range interaction. 
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The 2D scaling property of a finite range (Gaussian) in- 
teraction matches with maintaining fixed values of both 
the scattering length and the two-body bound state en- 
ergy. Thus, both mean- field zero-range and finite-range 
correlation calculations can use the strength determined 
to reproduce a given value of the scattering length in 
the Born approximmation. This is not possible in three 
dimensions, see the renormalization paragraph below. 

Radial solutions. To solve the hyperradial eq. (jSJ we 
need A. For p = only is present in eq. The free 

angular solutions are 

= P^o.A'-3) (^og 2a) , = + N-2), (24) 

where Pi"'^ '^^ is the Jacobi function Q and v — 
0,1,2,3... is a non-negative integer. For very small p 
perturbation then gives X^xi^'^ +2aBN{N ~l){p/hf. 

Increasing p for a repulsive or weakly attractive poten- 
tial leaves approximately the free wave function, but the 
energies are shifted from the contribution of the interac- 
tion. In first-order perturbation with the free wave func- 
tions the eigenvalues are denoted x'i\ where the lowest 
for p':^ b is, found to be 



A 



i/=0 



aBN(N - 1){N -2) 



(25) 



When the attraction is stronger, both wave function and 
energy change. However, the only p dependence in the 
angular equation is in v and in vi when a ~ 7r/2. Still, 
the eigenvalues A are essentially p independent when 
b <^ p <^ exp(— l/as), but perhaps lower than X,^. The 
eigenvalues must be obtained by numerical calculations. 

The behavior of the angular eigenvalues in the inter- 
esting range of p values are illustrated in fig.HJ They can 
be parametrized rather well by 



A(p)« 



xif^ 



ip/por 

2aBN{N - 1) 



iV- 2 ' 



(26) 
(27) 



where the last equality only holds for v = 0. Thus, po 
is independent of and inversely proportional to N. 
In these expressions we must for somewhat stronger at- 
traction replace A^f^ by the solution obtained numeri- 
cally. The curve for ob = —0.0559 in fig. |21 suddenly 
decreases dramatically when p/b 3000. This reflects 
the pronounced effect of the attractive pocket for large a 
appearing in fig. [21 for large p. 

The radial potential in eq. © has the same simple 
structure as the mean- field energy in eq. (|19|l . Analo- 
gously stable solutions only exist for constant A when 
A > -(TV - 3/2)(iV - 5/2). Using a[£o this implies col- 
lapse when OB < -{N - 3/2)(iV - 5/2)/[iV(iV - l)(iV - 
2)] « —1/{N + 1). This condition is again very simi- 
lar to the mean- field condition in eq. (|20|l . but now less 
surprising since the hyperradial and the mean-field radial 
equations both describe the same overall size dependence. 
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FIG. 2: The lowest angular eigenvalues as functions of p for 
different values of as and A^. The thick, solid line is Ao for 
= 20 and ub = —0.02. The parameters for the other curves 
are given directly on the figure. 



In fact a constant (j> in eq. H1U|I corresponding to = 
in eq. (|24|) leaves the total wave function as a function of 
only the hyperradius p. For Gaussian single-particle wave 
functions this is precisely identical to a Gaussian p depen- 
dence of /(p). For all other than Gaussian radial shapes 
of single-particle and hyperradial wave functions the con- 
nection cannot be made explicite. In the present formu- 
lation correlations, and thereby clear-cut deviations from 
the mean-field approximation, are included through the 
angular dependence of <&. These effects beyond the mean- 
field are perhaps most striking by the ability to describe 
simultaneously both states resembling a condensate and 
two-body bound states within the A'^-body system. This 
only requires two adiabatic potentials. 

For p values confined by the external field, we can in- 
sert the constant xlf ^ given by eq. H26I) in eq. Q which 
then is the radial harmonic-oscillator equation for an ef- 
fective angular-momentum quantum number /* defined 
by l*{l* + 1) = A + (iV - 3/2)(7V - 5/2) or equiva- 
lently I* = -0.5 + [A + (A^ - 3/2)(Af - 5/2) -h l/4]i/2. 
When /*(/* + 1) > -1/4, the solutions are well de- 
fined and characterized by the corresponding real val- 
ues of I* > —1/2. The eigenf unctions and eigen- 
values are (p/^O^'if"'^"^ [(p/^t)"] exp[-0.5(p/bt)2] and 
huj{2n + I* + 3/2), where L„ are generalized Laguerre 
polynomials and n is a non-negative integer. The ground- 
state energy is then fia;(/* -1-3/2) and the level spacing is 
precisely 2huj as for harmonic-oscillator excitations main- 
taining parity. 
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Renormalization of the interaction strength. In mean- 
field iD zero-range calculations the employed interac- 
tion is renormalized as ATTas5{r)fT? /m, where is the 
true s-wave scattering length derived from the low-energy 
two-body scattering properties. A short-range Gaussian 
representation of the delta function with a strength ad- 
justed to reproduce the same scattering length leads 
to totally different mean-field results. In order to get the 
desired results in the mean-field approximation with a 
Gaussian interaction, it is necessary to renormalize the 
Gaussian strength such that the Born approximation to 
the scattering length is reproduced instead of the true 
scattering length [13 ■ Thus the interaction in 3D mean- 
field calculations must reproduce the correct scattering 
length by using the Born approximation. 

This renormalization procedure guarantees three- 
dimensional mean-field results in agreement with the use 
of the above renormalization of the zero-range interac- 
tion. In the hyperspherical adiabatic method the large- 
distance asymptotic behavior of the angular eigenvalue 
should approach the expectation value for the mean- 
field zero-range interaction using the free wave function, 
i.e. the proper renormalization is seen to be equivalent to 
A,y for large p. Thus, to the extent the large-p be- 

havior of Aj/ is well described by aI'''' , we conjecture that 
this renormalization is not needed in 2D calculations. 

Let us now average a iD Gaussian interaction over the 
z-coordinate. Then we find that the Born approxima- 
tion to the scattering length of the new 2£'-interaction 

is as cx ag /hz, where 6z is the trap length in the z- 
direction and h <^ hz 0- Considering refs. 0, 0| 
and neglecting the logarithmic terms in their results, the 
present result is precisely the same for zero-range inter- 
actions when the true scattering length is substituted by 
the Born approximation. Thus, in ZD mean-field calcu- 
lations the Born approximation to the scattering length 
of the interaction potential, e.g. a Gaussian, is replaced 
by the true scattering length, whereas no replacement is 



needed in two dimensions. 

The stability condition, Nas > —1, is equivalent to 
Na'g /hz > —1, which coincides with the condition ob- 
tained by variation of the Gaussian length parameters for 
a deformed trap for a zero-range interaction where the 
3D scattering length is substituted by its Born approxi- 
mation 17]. Again this indicates that renormalization is 
not needed for 2I?-calculations. 

When the 3D scattering length is small compared to 
bz, the 2D and 31? coupling strengths for a zero-range 
interaction are related by 1720 ~ g3D/{bzV2TT). When 
the axial confinement is strong, i.e. bz «C a^, g2D becomes 
density dependent p^ . 

Treating correlations with the hyperspherical adiabatic 
method seems to require renormalization only in three 
dimensions where the mean-field and correlated results 
only coincide for differently renormalized interactions. 
For 2Z?-calculations the same interaction seems to pro- 
duce the same large-distance mean-field and correlated 
results. Thus, no renormalization seems to be required 
in 2D. 

Conclusions. We investigated symmetric iV-boson 
systems in the rigorous two-dimensional limit for attrac- 
tive two-body interactions. We derive the stability con- 
dition where the product of particle number, strength, 
mass, and square of the range has to be sufficiently small. 
The condition is the same for both mean-field and cor- 
related solutions and independent of the external field. 
Diatomic recombination into two-body bound states is 
not likely, since the two-body bound state would extend 
beyond the trap length for potentials sufficiently weak to 
allow stable solutions. Thus, these weakly attractive po- 
tentials, inevitably binding two particles, are able to sup- 
port stable condensates. These features are completely 
different from three-dimensional properties. These rigor- 
ous two-dimensional results cannot be compared to other 
calculations where the effective interactions are repulsive. 
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